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Abstract 

We present a dynamical alternative to inelastic dark matter as a way of reconciling the modu- 
lating signal seen at DAMA with null results at other direct detection experiments. The essential 
ingredient is a new form factor which introduces momentum dependence in the interaction of dark 
matter with nuclei. The role of the form factor is to suppress events at low momentum transfer. 
We find that a form factor approach is most likely not viable in the context of the standard halo 
model, however it is consistent with halo models suggested by recent Via Lactea simulations. As 
an example of possible form factors, we present a class of models where the necessary momentum 
dependence arises from interference of GeV mass gauge bosons coupling the dark matter to nuclei. 
At energies relevant for direct detection experiments these models contain one or two additional 
parameters beyond the case of a standard WIMP. 



1 Introduction 

For several years now, the Dark Matter (DAMA) collaboration has observed an annual modu- 
lation in the nuclear recoil event rate at its Nal detector, and the significance of the detection 
is now at the 8a level [1, 2]. It is striking that the modulation reaches its maximum within 
a week of June 2nd, when the earth's motion into the galactic dark matter is at its greatest. 
However, if one supposes that the DAMA signal is to be explained by a weakly interacting 
massive particle (WIMP) elastically scattering off of nuclei in the detector, then one pre- 
dicts far too many events to be seen at other experiments. The tension is particularly acute 
at low nuclear recoil energies, where standard WIMP event rates are expected to increase 
dramatically. 

Although DAMA apparently cannot be reconciled with other experiments within the con- 
text of minimal WIMP models (see, for example, [3]), this is not necessarily true in more 
general scenarios. Essentially, there is no other experiment that is sensitive to the DAMA 
signal in a model- independent context. In particular, there are three crucial aspects of the 
DAMA experiment which may be exploited in attempts at an explanation. These are: 1) the 
masses of the particular nuclei used in the experiments, 2) the fraction of the modulated vs. 
unmodulated event rate, and 3) the particular ranges of nuclear recoil energies being probed. 
Recently, inelastic Dark Matter (iDM) [4, 5, 6, 7, 8, 9] has received a great deal of attention 
for taking advantage of all three through a simple and elegant kinematic mechanism that 



easily arises in specific particle physics models. In fact, as the data has improved, few if any 
other models are consistent with DAMA and all null experimental results. In this paper, we 
will present an additional class of models as an alternative to iDM, which uses a dynamical 
rather than a kinematic mechanism. 

We will show how our mechanism can arise from a fairly straightforward model where the 
dark matter interacts with the Standard Model by containing constituents charged under new 
dark forces that kinetically mix with the photon. In terms of model-building, iDM is clearly 
simpler than our model. However, we wish to emphasize that model-building dynamical dark 
matter form factors can be an open direction for explaining DAMA, and there may be a 
variety of relatively simple alternative models. Thus, in the first part of this paper, we will 
try to determine in a model independant way how well one can possibly do with a form factor 
in terms of explaining all of the available data. We will see that the answer depends to an 
extent on the uncertainties currently present in our understanding of the dark matter halo 
velocity distribution. The picture works somewhat poorly in the context of the standard halo 
model, but much better with the results of the Via Lactea numerical simulations. In all cases 
we find a preferred dark matter mass window, lying in a general range from about 30 GeV 
up to about 50 GeV, with uncertainties coming from halo assumptions. 

The outline of this paper is as follows. In section 2, we will review the general ingredients 
which go into predictions for direct detection experiments. Section 3 will discuss the moti- 
vations for the form factor mechanism, as well as the model independent, best-case-scenarios 
for the picture. In section 4, we will present specific constructions which yield form factors 
with the appropriate momentum-dependence. We will note that the models we present sim- 
plify considerably in the low energy theory relevant for direct detection experiments. In this 
limit, the relevant parameters for the models become the mass of the dark matter particle, 
the overall cross-section, and the relative size of the q"^ term in the form factor compared to 
the term (as well as possibly higher-power terms). In this way, the simplest theories we 
consider have exactly as many free parameters as iDM. Section 5 contains a summary of all 
of the present experimental data relevant to our analyses, while in section 6 we will use this 
data to constrain the paramaters of our models, as well as give some examples of predicted 
spectra. Conclusions and future directions will be discussed in section 7. 



2 Review of Standard Elastic WIMP and Inelastic Dark 



Here we will briefly review the standard formalism for direct dark matter detection experi- 
ments; for detailed reviews, see e.g. [10, 11]. The rate per unit detector mass per unit recoil 
energy for dark matter nuclear scattering events is given by the general expression 
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Here Nt is the number of target nuclei per unit detector mass, and is the local halo 
density of the dark matter particle, whose mass is mDM- The rate depends on the local 
distribution f{v) of the dark matter velocity v relative to the earth. The integral over the 
velocity distribution depends on the nuclear mass and recoil energy only through the minimum 

velocity t'min = ^ \J "^^2'"' ^^^^ ^ dark matter particle must have in order to result in a nuclear 
recoil energy of Er, where mjv is the nuclear mass, and is the dark matter- nucleus reduced 
mass. 

The cross section for a standard WIMP on a nucleus with charge Z and atomic number 
A is given by 



da niN (Tp {fpZ + fn{A - Z)) ^ 



dEn 2v' pl /2 



F'n{Er). (2) 



where //„ is the dark matter-nucleon reduced mass, Up is the dark matter-proton cross section 
at zero momentum transfer, /„ and fp are the relative coupling strengths to protons and neu- 
trons, and Fn is the nuclear form factor. We will normalize our results to /„ = throughout, 
corresponding to the specific models we present. A common parameterization of the nuclear 
form factor is the Helm form factor 



where s = 1 fm, tq = \/r^ — 5s^, r = 1.2^4^/^ fm, and q = y/2rriNEji is the momentum 
transfer. A more accurate parameterization comes from the Woods-Saxon form factor, which 
is the Fourier transform of the Two-Parameter Fermi mass distribution p{r) — ^(r-'i)ia^y 
will take the Woods-Saxon form factor for Iodine (c = 5.593fm, a = 0.523fm), Germanium 
(c = 4.45fm, a = 0.592fm) and Tungsten (c = 6.51fm, a = 0.535fm) [12, 13], and use the Helm 
form factor for the remaining elements. 

Uncertainties in the form of the velocity distribution jiv) remain significant, as we will 
discuss in more detail, but a reasonable starting point is that /(w) in the galactic rest 
frame is a Boltzmann distribution oc e"^^/''^, with v = 220km/s, and with a cut-off at 
the escape velocity v^sc ~ 550km/s ± 50km/s [14]. In our analysis, we will take Vesc — 
600km/s, since this end of the allowed range tends to be more favorable for our models. 
To obtain our observed velocity distribution, the above expression must be boosted into 
the Earth's rest frame. The Earth is traveling through the galaxy with velocity Ve — Vs+ 
29.79km/s(cos(27r(t-tjuNE){0.262, 0.504, -0.823}+ sin(27ri-ijuNE){-0.960, 0.051, -0.275}), 
where the sun's velocity is Vs = {10, 225, 7}km/s, and we use coordinates with x pointing to- 
wards the center of the galaxy, y in the direction of the disc rotation, and z towards the 
galactic north pole [10, 7]. The magnitude of the Earth's velocity through the galaxy is at a 
maximum ont = tjuNE = June 2, causing a maximum in the net flux of dark matter particles. 

The inelastic dark matter scenario involves a simple generalization of the above such that 
a small amount of energy S, usually a mass splitting, is lost in the scattering to the dark 
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sector. The only modification of the above is to increase the minimum velocity by an amount 
5lq: 



^min - I {— + — ) + - (4) 

The additional contribution means that f min gets pushed out to the tail of the velocity dis- 
tribution for smaller momenta q and larger dark matter masses than in the elastic case. 
Its form is also very significant: ^ grows at decreasing q, and in fact Vy^im has a minimum 

possible value y^25~(mJ^^^~+^?T^^ that grows with decreasing nuclear mass. This is one of the 
main original motivations for inelastic dark matter. 

Event rates modulate the most over the course of the year for recoils dominated by particles 
on the tail of the velocity distribution. DAMA is the only experiment looking at a modulation 
signal, and ideally one way to try to reconcile DAMA with null results at other experiments 
would be to choose the dark matter mass sufficiently hght(e.g. [15, 16, 17, 18, 19]) in order to 
make Vmin at DAMA sufficiently large. Unfortunately, in the standard WIMP scenario, this 
means that the event rate at DAMA should rise rapidly as the recoil energy decreases [20] , and 
this does not match the observed spectrum. Furthermore, one faces immediate problems at 
the Germanium-based CDMS experiment, where the lighter element and lower recoil energies 
imply lower q and thus much less kinematic suppression. 

Fortunately, both of these problems stem from the sharp rise in the predicted spectrum 
at low q, and may be addressed by an additional form factor coming from the dark matter 
interactions, as we now discuss. 



3 Overlap in 

The tension between the apparent signal at DAMA and the direct detection experiments 
claiming null results becomes clearer through a comparison of their overlap in momentum 
transfers q, as shown in Fig. 1. In this paper, we will be considering scenarios in which 
the cross section eq. (2) is multiplied by an additional function -Fdm(q')^ coming from the 
dark matter sector. An important aspect of any such additional form factor -FbivKo') is that 
it does not depend on the nuclear mass and the recoil energy separately, but only on their 
product uinEji = q'^/2. As a consequence, any dynamical suppression at lighter nuclei can 
be completely compensated for by looking at higher energies. Perhaps more importantly, a 
form factor cannot be used to suppress events at a given experiment in an arbitrary way: the 
DAMA modulating signal probes values of q between qi ~ 70 MeV and qf ~ 120 MeV, and in 
this range one therefore has significant constraints on the nature of the form factor. Between 
qi and qf, the form factor cannot be too small lest the DAMA signature disappear, whereas 
outside of qi and qf, there is no model-independent lower bound. 

In practice, one would ideally like the form factor to serve two distinct, but related pur- 
poses, alluded to above. First, as noted, one would like it to suppress the scattering rate at 
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Figure 1: Overlap in q of the DAMA signal with several null experiments. The height of the 
null experiments has no particular meaning. 



momentum transfers outside of the DAMA range. In addition to this, it should also serve 
to flatten the spectrum of events observed at DAMA compared to the steep rise at low re- 
coil energies predicted by a standard light elastically scattering WIMP. Both of the above 
purposes may be served simultaneously by a form factor which falls off appropriately at low 
momentum transfers. 

Now, a key point is that the events seen by DAMA between and g/ essentially lead to a 
direct prediction (up to modulation fraction) for the events to be seen at other experiments 
within that same range of momentum transfers. These predictions are more or less indepen- 
dent of the choice of form factor, and it is therefore not immediately obvious whether they 
alone are enough to rule out form factor dark matter as an explanation for DAMA. The most 
basic question we must answer is thus the following: does there exist any function F{q) for 
the form factor - which we may take to be zero outside of the range Qi < q < Qf - which allows 
for the DAMA modulating signal, but which does not overpredict the number of events to be 
seen between and g/ at other experiments? 

Later we will consider explicit models that give rise to form factors, but for the moment 
we would like to answer this question while being as agnostic as possible about the model- 
building aspects. Thus we will begin by working with a physically unmotivated form factor, 
chosen solely with the goal of fitting the DAMA observed spectrum while simultaneously 
being consistent with the null experiments. To achieve this, we will construct a form factor 
to explicitly put the signal just below the la error bar at DAMA, bin-by-bin^. Furthermore, 
outside the range of the DAMA signal (i.e. below g = 80 MeV), we set the form factor 
to zero. An example is shown in Fig. 3. To evaluate the consistency of this form factor 
with experiment, we calculate the probability of the low number of potential signal events 

^More precisely, we construct the signal to be 80% of the signal-minus-lcr rate. 
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Figure 2: Schematic Feynman diagram for a form factor. 



at CDMS and CRESST-II using the pmax method [5, 21], which is based on the number of 
events predicted to have been seen in the gaps between potential events ^. We also calculate 
the of the fit at DAMA, which changes with dark matter mass and halo parameters due 
to the finite escape velocity. In particular, sufficiently large recoil energies necessarily have 
velocities above Wesc, and thus the signal is killed by phase space regardless of the size of 
the form factor. This tends to rule out small dark matter masses. We will take a fairly 
lax criterion for the probability of the DAMA fit in the lowest 12 bins, calculating only the 
goodness-of-fit for with 10 degrees of freedom. Of course, currently the form factor has a 
very large number of free parameters, but we are imagining the best-case scenario where one 
has constructed a model that automatically gives rise to the desired shape, and the only free 
parameters are the dark matter mass ttidu and the overall size of the cross-section (T„. 

The range of masses available depends highly on the choice of dark matter halo model, 
so let us now take a moment to describe our approach to this additional uncertainty. As 
we have mentioned, a commonly used simple model is a Maxwell-Boltzmann distribution for 
the dark matter velocities with a cut-off Vesc in the galactic rest frame. We will not attempt 
a systematic analysis of the effect of halo uncertainties, though in light of their effect we 
believe such an analysis would be very useful. Instead, we will follow [6] and consider two 
additional halo models, VL220 and VL270, based on the Via Lactea-I simulation [22]. Such 
an analysis is conservative in that halo parameters are notoriously difficult to constrain and 
a thorough analysis marginalizing over all halo parameters would open up more of the dark 
matter parameter space^. 

The VLa, models are based on a parameterization of the radial and tangential velocity 

^We present in detail our assumptions about the data from CDMS and CRESST-II in section 5. 

■^We emphasize that the main effect of the Via Lactea halo models will be due to their tighter velocity 
distributions, and similar qualitative effects would result from Maxwellian distributions with smaller RMS 
velocities. 
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Figure 3: An example of a spectrum from an idealized form factor. The discontinuous drop 
at low energies is by construction, to suppress all events outside of the DAMA signal range, 
while the discontinuous drop at high energies arises from the fact the escape velocity f esc cuts 
off the phase space beyond that point. 



distributions [3, 6] in the galactic rest frame, fit to Via Lactea simulation results: 

1 



fivR) = ^exp 
fM = ^r^exp 



2 \ 0'R-\ 

(5) 
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(6) 



with an = 1.09, ar = 0.73, vr = 0.72a/— t/ (ro), vt = 0.47a/— t/ (ro), where U{ro) is a normal- 
ization parameter. The resulting \/ —U (ro) around ro ~ 8 kpc in the Via Lactea-I simulation 
was 270 km/s; taking this value gives the VL270 model. If one instead sets \/ —Uir^ to the 
value 220 km/s in order to obtain a somewhat tighter distribution, then the result is the 
VL220 model. Given the uncertainty in the halo parameters and the difficulty in simulating 
the Milky Way (e.g. including baryons), we believe VL220 is within the allowed halo model 
parameter space and is important to include. 

The results of our idealized form factor construction are shown for these three different 
halo models in Fig. 4. One can see that, even with our fairly conservative constraints, there 
is hardly any region of masses allowed within the standard halo model. The rapid rises in 

for the fit to DAMA occur when the minimum velocity in an energy bin starts to exceed 
the bound from the escape velocity, at which point no form factor is capable of reproducing 
the DAMA signal in that bin. This suggests that a purely dynamical explanation of DAMA 
is disfavored in the standard halo model^. However, there is significantly more room to work 



''The range of allowed masses is somewhat sensitive to the value of the quenching factor q\ for iodine 
recoils at DAMA, as it affects the reconstructed nuclear recoil energy range of DAMA's signal. Here we 
have taken q\ = 0.085, but if future improved measurements find it to be lower, then that would move 
the DAMA spectrum up in recoil energy, thereby weakening the CDMS constraint and strengthening the 
CRESST constraint. For example, in the standard halo model, qi — 0.08 opens up the 95% allowed mass 
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Figure 4: Constraints on idealized form factors for different halo models. For each dark matter 
mass, we choose a form factor to explicitly match the DAMA spectrum, as described in the 
text. The confidence levels at DAMA, CDMS and CRESST2 are then plotted. We use a 
goodness-of-fit test to define the confidence at DAMA, and the pmax method for CDMS and 
CRESST2. The standard halo model leaves very little room for a form factor explanation for 
all of the data. Highlighted regions have all constraints below 95%. 



with in the VL270 model and, especially, the VL220 niodel^. 

At this point, it is worth re-emphasizing some of the main differences between the iDM 
mechanism and a dark matter form factor. Because the effect of the splitting 6 in eq. (4) 
grows at smaller q, one obtains a very rapid kinematic shut-down at low DAMA energies, 
roughly of the form ~ e~^^^''^^ . One of the difficulties of model-building a form factor will be 
to obtain a sufficiently rapid shut-down purely from the dynamics. A more drastic difference 
is that fmin is a monotonically increasing function of q for elastic dark matter but not for 
inelastic. This makes it generally impossible given the recoil energies and nuclear masses 

range to 30.5GeV < moM < 35.5 GeV. However, at qi = 0.07, the CRESST constraint is stronger than CDMS 
and the mass range shrinks to SBGeV < ttt-dm ^ 35 GeV. 

^We have checked that scattering off of sodium does not work in any of the three halo models we consider 
for any range of dark matter masses. 
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at DAMA and CDMS to choose parameters so that there are no events at CDMS while 
simultaneously maintaining the DAMA signal. However, as has been noted in the literature, 
this is a possibility for inelastic dark matter if the splitting 5 is chosen sufficiently carefully. 
We can sharpen this prediction somewhat by considering the minimum number of events 
predicted at CDMS by our idealized form factor in the 95% goodness-of-fit DAMA region 
(i.e. Xdama < 18.3). This number is proportional to the effective exposure ecDMS at future 
CDMS runs, so that 

AT eCDMS 

,min '"CDMS.min f^^^i -t i \' ) 

200kg day 

where we find that rcDMS,min ~ 4.4,1.1, and 2.2 for the standard halo, VL220, and VL270 
models, respectively. 



4 Models 

Here we present models leading to a form factor which is consistent with present data. One 
very simple possible form factor, which is relatively generic at sufficiently low energies, is 
F{q) (X (f. Such an interaction results in a factor of only in the event rate, and leads 
to insufficient suppression of events at null experiments. We are thus lead to introduce a 
mechanism to further reduce events outside the DAMA region. 

Our basic idea is to have multiple gauge bosons coupling regular matter to dark matter 
which then naturally interfere with one another. In particular, we will examine two models. 
The first describes two gauge bosons and leads to an effective form factor 

The second model includes three gauge bosons, leading to 

W)~(^-2^ + ^). (9) 

In both cases, we will see that the models possess natural regions in parameter space 
where there is a custodial symmetry relating the gauge boson masses. In those regions, the 
analysis of the models simplifies. Specifically, choosing all mj to be near a GeV, at momenta 
of interest {q <^ rrii), 

F2GB{q') = C2q\q^-ql) (10) 
F,GB{q') = m\q'-qlW-ql). 

Here, the parameters qi are of order 10 MeV, while the q are absorbed into the overall 
crossection. Obtaining a crossection of the appropriate size in both cases requires A, ~ 
GeV. Therefore, as far as the direct detection experiments are concerned, the two and three 
gauge boson models contain two and three parameters, respectively. 
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4.1 A two gauge boson model 

(1 2) 

Consider two massive gauge bosons Aj^ ' which mix with hypercharge with opposite signs. 

^-e{grF!,])-g,F('J)B^'^ (11) 

Such a mixing can arise, for example, from integrating out two heavy fields each of which 
couples to one of the gauge bosons and to hypercharge with charges chosen appropriately. 

We assume that the dark matter particle is not charged under the gauge bosons, but that it 
may be made of constituents that are charged. Consequently, it develops higher-dimensional 
interactions with both gauge bosons. If it is a scalar, the leading higher-dimensional interac- 
tions take the form 

^ - ^ {jt^F(]) + f^F(^j) d,X*d.X. (12) 

Here X is the dark matter particle and A is characteristic of its size. Thus, the exchange of 
the two gauge bosons between the dark matter particle and the protons in nuclei will lead to 
the form factor of Eq.(8). 

Let us now discuss how the gauge bosons obtain their masses. As mentioned, the form 
factor simplifies in the limit that a custodial symmetry relating the gauge boson masses is 
weakly broken. We will therefore describe the Higgs fields and their potential with that in 
mind. Consider a complex doublet, $ = (0i,02)- The custodial symmetry is an 5*0(2) 
rotation of 0i into 02- ^ couples to the gauge bosons as 

i^,$=(a,-^|A«-^|4^)a3)$. (13) 

We now turn to the Higgs potential. Let us assume in the following that the only source of 

(2) 

breaking of the SU{2) symmetry acting on $ comes from interactions with A/j , and that 
the dark sector interactions respect a charge conjugation symmetry (exchanging 0i and ^2) 
which takes A^^^ — >■ —A^^\ With these symmetries, the $ potential takes the form 

y($) = -m2$^$ + A($^$)2 + y($V3$)^ (14) 

If y > 1, at the minimum $^3$ = 0. This ensures that the gauge boson mass matrix is 
diagonal and that there is an SO{2) custodial symmetry relating the gauge boson masses, 
such that mi^2 = 91,2 v. Including terms in the potential which violate the charge conjugation 
symmetry would make this potential generic, resulting in general gauge boson masses. For 
our purposes, however, it is useful to assume the there is only a small violation of the above 

2 2 

custodial symmetry, such that ^ — = (dm)^. Such a violation could come, for instance, 
from the above mixing of the gauge bosons with hypercharge (which violates the charge 
conjugation symmetry at order e^). In this case, the size of (5m)^ would depend on additional 
physics such as the cutoff scale of the $ theory, ((5m) ^ ~ jJ^A^^^^jj. One can imagine other 
sources for {Sm)^ as well, and we therefore leave it as a free parameter. Thus, in the limit of 
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interest, where (f <^ Tn\ 2-i the form factor is approximated by 

F,GB{q') - ^ f ^ - {5mA ^ cn\i - ql). (15) 

" \ yi92 J 

We note that besides the constant c that can be reabsorbed into the total crossection, the 
only relevant parameter for dark matter detection experiments is gfo- 

4.2 A three gauge boson model 

One can also consider extending the above scenario to include an additional gauge boson which 
mixes with hypercharge. The interference between the three gauge bosons can naturally 
achieve further suppression of the cross-section at low momentum transfer. The resulting 
model is similar to the model above, but with an 5*0(3) custodial symmetry instead. As 
before, the three gauge bosons mix with hypercharge, 

L ^ e (^iFt;,) - Ig^Ff} + ^3^^'^) B>^r (16) 

Such a mixing can be obtained by integrating out three heavy fields, each charged under 
hypercharge and under only one of gauge bosons. All heavy fields have the same charge 
under their respective gauge bosons, and charges +1, —2, and +1 under hypercharge. 

The dark matter particle couples to the gauge bosons through the leading higher-dimensional 
interactions, 

£ = ' d'-^i' + %P^S + fl^-^?) 9.A-a„x (17) 

Hence, gauge boson exchange results in the effective form factor of eq.(9). 

We will again concentrate on a version of this model with approximate custodial symmetry. 
To that end, we take the Higgs field, $, to be a four of S\J (4), charged under the gauge bosons 
as 

D^^ = {d^ - iQiAfQ^ (18) 
Here, the charges, Qi are embedded in SU (4) as 

where 1 is the two by two identity matrix. Assuming that the only interactions breaking 
the SU (4) symmetry are the gauge interactions, a general potential respecting all manner of 
charge conjugation symmetries is 

V($) = -m2$t$ + |/,($tg^$)2 + ^($tg2g3$)2_ (20) 

If > 0, then at the minimum an SO{^ custodial symmetry acting on the gauge boson 
masses will be respected, and = giv. This custodial symmetry will be violated at order e^. 
We thus define ibrn)^ — — ^ — 
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In order to obtain a simplified version of this model at low q, we will further require 
that the three gauge bosons unify at some scale Mu. For instance, the unified group can be 
SO{S) D SU{4:) X C/(l), with $ transforming as the eight of SO{8). Therefore, gi{Mu) = gu- 
Below Mu, we assume that the matter content is such that the running splits the couplings, 
but maintains the relation ^ — ^ + ^ = 0. This is possible if for example there are twice as 
many fermions charged under only as under only U{1)2, with no new matter (besides 

$) charged under f/(l)3. Threshold effects would then account for a violation of the relation 
between couplings, so that — + ^ = Sg < jj^^. Consequently, at q <ti rrii, 

9l 92 93 ^"-"J 



1,6 



4 4 ' 4 

91 92 9l 



+ \ ]q'^ -Sg v'^q^ + {Smfv'^ 



Csq'iq'-qW-qi). (21) 



We note that choosing Acutofr ~ 477^; naturally gives ((5m)^v^ ~ 10 ^v^, which is the appro- 
priate size for qf (assuming all couplings are of order one). 



5 Null Experiments and Analysis 

We will briefly describe the null experiments whose constraints we include and the assumptions 
that go into our analysis. We have attempted to follow the analysis assumptions of [5] in 
order to make the comparison with inelastic dark matter transparent. Summary tables of the 
assumptions at each experiment are 1,2, and 3. Tables 2 and 3 show the relevant data used in 
our analysis concerning the various null experiments. Generally, these experiments see some 
set of unexplained events whose energies are indicated in table 3, and in our analysis we have 
allowed ourselves to interpret these as potential dark-matter-induced nuclear recoils. The 
energies listed in this table are actual nuclear recoil energies, and have thus been obtained by 
dividing the detected energy by a quenching factor where appropriate. In all null experiments 
in which individual potential events are published, we calculate the constraints using the 
pmax, or "maximum gap", value [21]. This is the probability, given a specific model, of 
an experiment having seen no events in the gap between any two adjacent potential events. 
When the number of observed events is very low, this method is very sensitive to, for example, 
non-blind cuts on events, since a single additional event in the middle of a gap can drastically 
increase the hkelihood. 



5.1 DAM A 

We take for the DAMA signal the published modulation amplitudes from Fig. 9 of [1], which 
we summarize in Table 1. The energy shown is that observed in the detector, and is smaller 
than the actual nuclear recoil energy by an amount given by the "quenching factor". For 
Iodine scattering events in a Nal crystal, we take the quenching factor to be 0.085. We 
implement scattering off Iodine only, as our form factor suppresses events scattering off of 
Sodium. We calculate the modulation amplitude as half of the difference between the event 
rate on June 2 and Dec. 2. To obtain our confidence limits on our models from DAMA, we 
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nil _L nnoo 


6.5 ± 1. 


.074 ± .054 


f; o f; _i_ o f; 
O.ZO ± .ZO 


.UUoo ± .UUoo 


7.5 ± 1. 


-.0008 ± .037 


'tf; _i_ of; 
.ZO 


.UUzo ± .UUoo 


8.5 ± 1. 


.051 ± .036 


6.25 ± .25 


.0055 ± .0035 


9.5 ± 1. 


.051 ± .027 


6.75 ± .25 


.0013 ± .0033 


10.5 ±1. 


.065 ± .026 


7.25 ± .25 


.00025 ± .0035 




7.75 ± .25 


.00025 ± .0035 



Table 1: Left: amplitude of the modulating part of the DAMA signal, divided into 0.5keVee 
bins. Right: amplitude of the measured rate at KIMS, in 1.0 keVee bins. 

find the best at a fixed mass and show contours of = Xmin + 4.61(9.21) for 90% (99%) 
confidence. 

5.2 CDMS 

We combine the CDMS data from the three runs in [23, 24, 25], for a total of 174.4 kg days 
after including their pubhshed 30% acceptance rates on a standard WIMP. The first run had 
20 kg-days and ran from Oct. 11 - Jan 11, the second had 34 kg-days and ran from Mar. 
2 5- Aug 8, and the latest five-tower run (which saw no events) had 121 kg-days and ran from 
Oct. to July. 

5.3 CRESST-II 

We combine the CRESST-II data from the 2004 run [26] and the latest commissioning run 
[27] from 2007. The '04 run ran from Jan. 31 - Mar. 23, had 20.5 kg-days of exposure, and 
looked at recoil energies of 12-50 keV. The '07 run ran from Mar. 27 - July 23, had 47.9 
kg-days of exposure, and looked at 12-100 keV. Both runs had very mild cuts leading to a 
90% acceptance rate. The '04 run had two detectors, one of which (JULIA) saw four events, 
and the other of which (DAISY) saw no events using the published quenching factors, and 
one events after varying the quenching factors within their published uncertainties. The '07 
run also had two detectors, which saw a total of seven events. We should note that the pmax 
constraints from the '07 run alone are much stronger than the constraints from the combined 
'04/'07 analysis, but given that the '07 run was non-blind we follow [5, 28] and take all the 



13 





n .1 OTYl Oiri" 






Olgjlldi VVlllLlUW 

{Kev ) 






Uct. z- July oi 


1 '7 A A 


in 1 nn 






Uct. i - reb. Z5 


1 00* 


z - /U 


PRpqcT TT 'DA 




Jfill. 01 - iVlclI. ZO 


iO.O 
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CRESST-II, '07 


CaW04 


Mar. 27 - July 23 


43.1 


12 - 100 


ZEPLIN-II 


Xe 


May 1 - June 30 


112.5 


13.9 - 55.6 


ZEPLIN-III 


Xe 


Feb. 27 - May 20 


126.7 


10.7 - 30.2 


KIMS 


Csl 


see text 


N/A 


20.6 - 65.4 



Table 2: Elements, dates, effective exposures and signal windows for the null experiments. 
The exposures have been listed here only after efficiencies and cuts have been talcen into 
account (it is still necessary, however, to multiply by the fraction of target nuclei in the 
detector material). *XENON effective exposure is approximate, since the efficiencies and 
acceptance rates arc energy-dependent. This was talcen into account in our analysis. 

CRESST-II data together. 

5.4 XENON 

For XENONIO [29], we tak;e the published data of the recent reanalysis [30] in the extended 
range 2.0-70 kcV. The experiment ran from Oct. 1- Feb 28. and had a total of 316 kg-days. 
We take the published energy- dependent acceptance rates and efficiencies. 

5.5 ZEPLIN-II and ZEPLIN-III 

ZEPLIN-II is a Xenon based experiment with 225 kg-days of effective exposure after cuts are 
applied, and in addition there is a recoil detection efficiency of 50%. We take the quenching 
factor to be 0.36, and the 29 events from Fig. 8 of [31]. We take all events as potential signal 
in calculating our limits. 

ZEPLIN-III ran in 2008 from Feb. 27 - May 20 and observed only 7 events between 
10.7-30.2 keV [32]. 

5.6 KIMS 

KIMS [33] had four CsI(Tl) crystals looking in the range 3-11 keVee, two of which ran roughly 
year-round and two of which ran between June and March. We take the weighted average 
over all four crystals to obtain the measured rate and error, shown in table 1. We use the 
parameterization of the quenching factor from [6]: 

175e-^«/i37 

^^^^^^^^ ^ l + 0.00091i?«' ^^^^ 
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J-JVCiltto I l^C V j 


CDMS 


10.5, 64 


XENON 


15.8. 23.8. 21.L 25.9. 32.2. 33.7. 38.7. 39.2. 10.1. 13.3. 13.6. 51.6. 62.3 


CRESST-II '04 


13, 18.5, 22.5, 24, 33.5 


CRESST-II '07 


16.7, 18.03, 33.09, 43, 43, 43, 63 


r~71 — T1VT XX 

ZEPLIN-II 


13.9, 15.3,16.7,16.7,16.7.18.1, 
19.5,19.5,20.9,20.9,22.2,22.2, 
23.6,27.8,27.8, 29.2,30.6, 33.4, 
37.5, 37.5, 38.9, 41.7, 43.1, 44.5, 
44.5, 48.7, 50.0, 50.0, 54.2 


ZEPLIN-III 


14.6, 17.2, 20.8, 22.9, 23.4, 25.1, 28.1 



Table 3: Unexplained events observed by the null experiments, listed by nuclear recoil energy 
in keV. 



and define the constraint as the requirement that the total scattering off Cesium and Iodine 
not be more in any of the first five bins (3-8 keVee) than the measured rate plus 1.64 times 
the error. 

6 Results 

Exclusion plots are shown in Figs. 5 and 6^ for the two- and three-gauge boson (2GB and 3GB) 
models, respectively. In order to reduce the parameter space and thus simplify the contour 
plots, we have taken both nodes in the 3GB model to be at qi = q2 = qo; the improvement 
in consistency with the null experiments by varying them separately is relatively small. We 
have checked that in each case, the constraint from KIMS can be neglected since it is weaker 
than both the CDMS and CRESST-II constraints, and in fact its 95% constraint curve does 
not cut out any of the DAMA 90% preferred region^. The DAM A inner (outer) contours 
are shown at 90% (99%), and the pmax constraints for the null experiments are shown at 
99% (95%) for the 2GB (3GB) model®. There is a fairly strong dependence of the hmits on 
the dark matter mass, both at DAMA and at the null experiments. As discussed in section 
3, sufficiently low masses put the DAMA signal above v^sc and thus are excluded. On the 
other hand, sufficiently high masses cause the velocity dispersion to flatten out at low recoil 

^The dark matter-proton cross-section (7p now involves the form factor, and so is a function of q. We 
therefore have defined o-p^iooMeV in the foUowing way: ap = ctj,^ iooMeV-pDM(<?)> where FoMiq) is normaUzed 
e.g. as Q2(g2 _ ^2)(^2 _ g2^/(iooMeV)6. 

^We note that in our models the constraint from comparing the unmodulated DAMA signal to the modu- 
lated DAMA signal is an even weaker constraint on the modulation fraction than KIMS is, and so may also be 
neglected. The danger of overpredicting the unmodulated signal in low energy bins En < 20 keV is avoided 
as we predict the spectrum to decrease at low energies. 

^Thc reason for this difference is that the DAMA 99% preferred region for the 2GB model is just barely 
excluded at 95% by our analysis. We hope this will not cause too much confusion. 
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Figure 5: Exclusion plots for the two-gauge-boson (2GB) model. Confidence limits from 
DAMA are shown in purple (light blue) for 90% (99%), all other confidence limits are 99%. 
The spectrum at DAMA is shown for a benchmark point with VL220; see table 4 for the model 
parameters, is taken in units of cm^. 



energies. In our specific models, the lower energy bins are too suppressed by the form factor 
and the spectrum is incorrectly distorted. On top of this, larger masses do not lead to a large 
enough modulation fraction at DAMA in order to evade the null constraints. In Figs. 5 and 
6, we have chosen masses from around the regions where the constraints are weakest. 
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Figure 6: Exclusion plots for the three-gauge-boson model (3GB) with both nodes taken at 
the same value qi = q2 = qo for simplicity. Confidence limits from DAMA are shown in purple 
(light blue) for 90% (99%), all other confidence limits are 95%. The spectrum at DAMA is 
shown for two benchmark points; see table 4 for the model parameters. 
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The DAMA preferred region for the 2GB model is excluded by our analysis at 95% con- 
fidence. In all cases, the culprit is CDMS, due to the fact that the form factor does not fall 
quite rapidly enough below the DAMA signal region. We find though that at VL220 the 95% 
CDMS line just barely excludes the DAMA preferred region, and there are still a mimbcr 
of uncertainties (e.g. in quenching factors, energy resolution) that we have not attempted 
to take into account. The situation is much better in the 3GB model, where we show 95% 
exclusion limits. In the standard halo model, no parameter space is left at 90%, but for the 
Via Lactea models the DAMA preferred region is just barely allowed by our exclusion limits. 
The GDMS constraint is considerably weaker for the 3GB model, due to the sharper fall-off 
of the form factor at energies just below the DAMA signal. As one might have expected, 
the other experiments with a sizeable range below the DAMA signal (XENON, ZEPLIN-II, 
and ZEPLIN-III) have also gotten weaker relative to CRESST-II. One can clearly see the 
effect of the node qq as it passes through the bulk of each individual experiment, removing 
the predicted events in that region. For go above about 70 MeV (60 MeV) in the 2GB (3GB) 
model, however, the node starts to suppress the DAMA signal, whereas ideally the form factor 
will be more or less set in this region by the DAMA spectrum. Around ~ 80 MeV, the 
CRESST-II signal is drastically suppressed by the node. Thus the trend in both models is 
for the lower-Q' experiments to dominate the constraints at smaller q^, whereas for 40 MeV 
i$ ?o 70 MeV, they are suppressed and CRESST-II becomes the dominant constraint. 

The halo model with the most open parameter space is in all cases VL22o- The first 
major effect of the smaller v at VL220 is to tighten the velocity distribution. This allows for 
larger dark matter masses without sacrificing the necessary modulation fraction. Additionally, 
larger masses effectively push the escape velocity cut-off to higher recoil energies at DAMA, 
improving the fit to the upper energy bins. Finally, it is important to note that despite 
the tighter distribution at higher energies, the distribution in the lab frame is still relatively 
fiat below the earth's velocity, as the velocity of such events are dominated by the earth's 
contribution rather than the halo fiuctuations. Given how much a modest change (18.5%) in 
V can affect the allowed parameter space, it is clearly an important issue to understand the 
experimentally and theoretically allowed range of parameter space for halo models, so that 
it would be possible to more systematically include all such uncertainties in setting exclusion 
limits. In the absence of such an analysis, the precise meaning of exclusion limits is not 
entirely clear. 

The value of for the DAMA 90% and 99% contours depends on the minimum 
for each plot, so it can be helpful to know this value. For the 2GB (3GB) model, — 
4.6,4.7,4.7(9.8,5.0,5.24) for the standard halo, VL220, and VL270 plot, respectively. In all 
cases, x^in is below 10 = 12 — 2, as appropriate for 12 DAMA bins and 2 free parameters (up 
and qo). Given the strong dependence of the fits on moM, it is possible that the x^ distribution 
for our models more closely follows a distribution with three free parameters, in which case 
our Xmin is still sufficiently small at VL220, VL270, but the DAMA contours will be larger. In 
this sense also, our analysis is conservative. 

In table 4, we show three benchmark models, one for the 2GB model with the VL220 and 
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VL220 


2GB (go = 50MeV) 


50 GeV 


11.8 


2.34 X 10-^°cm2 


0.97 


0.89 


VL220 


3GB (gi = 42.5, ga = 38MeV) 


36 GeV 


8.9 


2.00 X lO^^^cm^ 


0.90 


0.90 


VL270 


3GB (gi = 50, ga = 37.5MeV) 


32 GeV 


10.3 


2.10 X lO-^^cm^ 


0.94 


0.95 




)7 


V 


14.9 


1.79 X lO-^^cm^ 


0.90 


0.90 



Table 4: Our three benchmark points. We also indicate what the constraints on the 
3GB(VL27o) benchmark point would be if a slightly smaller overall cross-section CTp^iooMeV 
were chosen. 




£„(kcV) £«(kcV) E„(keV) 



Figure 7: Predicted spectra at Ge (left), W (center) and Xe (right) for the three bench- 
mark points in table 4. We have chosen the time of year to correspond to the dates of 
CDMS, CRESST-II '07, and ZEPLIN-III, respectively. The 2GB(VL22o), 3GB(VL22o), and 
3GB(VL27o) models are shown in sohd, dashed, and dot-dashed. 



two for the 3GB, one for each VL^; model, along with their statistical constraints from CDMS 
and CRESST-II. Note that the null constraints are within 97% for the 2GB model and within 
95%(90%) for the 3GB, VL27o(VL22o) niodel. We have also shown in table 4 the constraints 
on the 3GB(VL27o) model using a slightly larger value for the cross-section CTp^iooMeV, because 
such a cross-section would be permitted by a goodness-of-fit test. More precisely, for 
the fit to DAMA in the new 3GB(VL27o) model is 14.9, which corresponds to a probability 
of 86.4% for 10 degrees of freedom, and 90.6% for 9 degrees of freedom. It is, however, just 
outside the 99% contour = xLin + 9-21 = 14.4, which underscores that the constraints 
are sensitive to the statistical assumptions being made. Finally, we also show in Fig. 7 the 
predicted spectra for our three benchmark points at CDMS and CRESST-II. 

7 Discussion and Future Directions 

In this paper, we have discussed the possibility of using a dynamical form factor to recon- 
cile the DAMA annual modulation signature with the null results of other direct detection 
experiments. We have shown that, independent of specific model assumptions, such a form 
factor can in fact be consistent with all available data. The picture works to varying degrees 
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depending on the assumptions one makes concerning the distribution of velocities in the dark 
matter Halo. 

An important question is how one might distinguish this scenario from that of inelastic 
dark matter. For this purpose, probably the most important indicator would be in the energy 
spectra to be seen at future runs of the various direct detection experiments. 

As discussed in section 3, independent of model details, form factor dark matter leads to 
a prediction for the events to be seen at any direct detection experiment within the range 
of momenta probed by DAMA. Although this prediction depends somewhat on the dark 
matter mass, some non-trivial event rates will always be present. Inelastic dark matter, on 
the other hand, has a bit more flexibility in this regard: due to the form eq. (4) of the 
minimum velocity for inelastic scattering, it is possible to choose parameters such that Vj^in 
is always above the escape velocity at CDMS but below the escape velocity at DAMA in the 
energy region of the modulation signal [4, 5]. The result is that inelastic dark matter may be 
consistent with possibilities such as zero events at future runs of CDMS, a situation which 
would rule out a pure form factor. 

In contrast to this, there are certain aspects of the form factor setup which are actually 
more flexible than in the inelastic case: inelastic dark matter makes a flrm prediction that all 
experiments should find zero events at momenta smaller than a particular cutoff gmin ~ ^ , 
which must be set approximately to the lower end of the DAMA spectrum. A general form 
factor, on the other hand, could accomodate such events if they were to be seen. 

Further differences between the predictions of inelastic and form factor dark matter could 
become relevant after large amounts of data have been collected. In particular, details of 
the energy dependence of the fraction of modulating events could be used to differentiate 
the two pictures, but only with large amounts of statistics. It is also possible that future 
directional detection experiments could provide an additional handle [8]. Finally, although 
collider signatures for the general form factor scenario are model dependent, the specific setups 
discussed in section 4 could lead to events at low energy e'^e" colliders of the types discussed 
in [34, 35]. 

There are a wide variety of future directions to pursue within the general framework we 
have discussed in this paper. Of course, we hope that future simulations of the galactic 
halo, including baryons, will yield a more precise picture for the structure of the dark matter 
velocity distribution. At the present level of understanding, uncertainties in this distribution 
have a significant impact on the possible viability and allowed parameter ranges for not only 
form factor dark matter, but for essentially any proposed explanation for the DAMA annual 
modulation signature. 

From the particle physics side, it would certainly be interesting to try to write down 
additional types of models leading to form factors with the general features discussed in 
section 3. In this paper we considered a simple "proof of principle", but other model building 
directions should be viable as well. Within the context of the specific constructions presented 
here, it would be worthwhile to try to flesh out a more detailed structure for the composite 
dark matter particles, along with a corresponding cosmological history. 
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Finally, considering form factor dark matter in the context of a possible "channeling" 
effect [36, 37] in the DAMA/Nal detector seems to be a very promising direction, and will be 
the subject of a future paper [38] . 

Note added: After this paper appeared on arXiv, there appeared [39], which considers 
certain dark matter form factors in the context of spin-dependent cross-sections. 
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